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The properties of a nonrelativistic charged particle in two-dimensions 
in the presence of an arbitrary number of nonquantized magnetic 



fSJ ' fluxes are investigated in free space as well as in a uniform magnetic 



field. The fluxes are represented mathematically as branch points in 
one of the complex® coordinates. To construct many-flux solutions 
without a magnetic field, however, the fluxes must be treated as dy- 
namical objects dual to the charges. A medium made up of fluxes acts 
like an anti-magnetic field and tends to expel the charges. 



Introduction 



^ ■ In a classic paper, Aharonov and Bohm have pointed out that a locally 

. trivial vector potential of a gauge field can lead to observable effects in a 

nontrivial topology through the phase of the wave function of a charged 
particle. This has spawned many important theoretical concepts in later 
developments in gauge theories, like instanton, monopole, and confinement 
by monopole condensation. On the experimental side, confirmation of the 
original Aharonov Bohm (A-B) effect due to a nonquantized magnetic flux 
tube has been somewhat controversial because of the difficulty of setting up 
ideal conditions, but there have been several experiments, notably those by 
Tonomura et al showing the shifts in diffraction fringes due to magnetic 
fluxes, as expected by the A-B effect. (See for a general review of the 
problem.) The present work was motivated by the theoretical question: What 
would be the properties of a medium filled with many nonquantized magnetic 
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fluxes (vortices)? Assume, for simplicity, that the fluxes are all parallel so 
that problem can be reduced to that of an electrically charged particle in a 
2-plane pierced by pointlike magnetic fluxes. Then, for example, what would 
be the behavior of a charged particle in a lattice of such fluxes? What would 
be its energy spectrum? Would the fluxes cause more drastic effects than 
a phase shift in the wave function? Although the semiclassical arguments 
suffice to discuss and analyze most situations realistically, it requires a more 
rigorous treatment of the Schroedinger equation to answer these questions. 
This turns out to be a highly nontrivial problem. To our knowledge, the 
solution of the Schroedinger equation in the presence of two or more fluxes 
has not been explicitly constructed. (A qualitative discussion of special cases 
was made by Peshkin et al. [Q. Some other theoretical problems related to 
the present one have been discussed by Aharonov and Casher 0, Dubrovin 
and NovikovH, Jackiw||ll||, and Lewisl 



We develop here a method of solving the many-flux problem by adapting 
the work of Sommerfeld|P [H. He solved the problem of diffraction of light 
by a semi-inflnite wall in two dimensions by regarding the wall as a branch 
cut in a two-sheeted Riemann surface. The incoming plane wave in the 
physical plane dives into the unphysical second sheet as it hits the cut, while 
the reflected wave emerges from the second sheet through the cut. The 
solutions were constructed as a contour integral of a kernel in a complex 
angular variable over the two sheets. 

The A-B problem can be posed in a similar fashion. The pure gauge vector 
potential of a flux is singular at the flux site. In the complex variables z = 
x^iyW behaves like {a/2){i/z, —i/z)^ where 27ra is the flux strength. After 
removing the potential by a gauge transformation, we get a wave function 
that satisfles a free Schroedinger equation, but it is singular in the sense that 
the function acquires a phase rj = exp(2z7ra;) after encircling the flux once, 
so the wave function must vanish at the flux if a is noninteger. Regarded 
as a function in the complex 2-plane, the problem then reduces to that of 
solving the free Schroedinger equation on a Riemann surface where the flux 
is a branch point, and the solution must satisfy the boundary condition that 
it gets a phase change t] when going around it once, and vanishes at the 
branch point. 

The paper is organized as follows. We flrst establish the basic formalism 
of complex integral representation, apply it to rederive the known results for 
the case of a single flux with or without the presence of a uniform magnetic 



fleld.[n|14|. Then the problem of many fluxes is addressed by two methods 
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which yield different types of solutions. One of them is to use multiple integral 
representations, and is applicable only if a magnetic field is present. The 
other is to treat the fluxes as dynamical objects having properties dual to the 
charges. An important general property that emerges is that the nontrivial 
phase information around all fluxes must be represented as branch points in 
either the z or the z coordinate, but not a mixture of them. It means that, 
since the wave function must vanish at the flux sites, the (fractional) angular 
momentum around each flux, as opposed to the intrinsic flux strength, has 
the same sign, either positive or negative. Thus for example, if the two 
flux strengths are ai > and 02 < 0, the solution is constructed as having 
positive angular momenta ai + /i > and /2 + «2 > 0, /j = integer, in the 
z variable (and similarly negative angular momenta in z). Otherwise the 
solution would have to be represented by a function in 2; in a region around 
one flux, and by a function in z in another, but it is not possible to match 
the values of the wave functions and their derivatives at the interface. Q 

The fact that angular momenta around all the fluxes must have the same 
sign is not intuitively clear, but may be inferred from the zero energy solu- 
tions, where the nontrivial information about all the fluxes must be encoded 
in terms of either analytic or antianlytic functions. It is also understandable 
from the following consideration. Suppose a\ = — 02, and let them approach 
each other and collapse to nothing. The wave function originally vanished at 
both sites and was finite elesehwere. So it will vanish at the point of collapse, 
hence it will be in a p or higher angular momentum state. 

One of the consequences of the above property is that a medium made 
up of fluxes tends to expel the charge as if it was placed in an anti-magnetic 
field. This behavior has a simple explanation. By construction the angular 
momentum around each flux is of the same sign, say a', irrespective of the 
sign of its intrinsic flux strength. The total angular momentum around a 
circle at a radius R lattice units will be ~ a'lrR'^, hence the radial momentum 
p{R) is ±{E — (a'TrRYY^"^ . For large R, then, the wave function will go like 

^To see this, consider at zero energy two fluxes of opposite sign located at ±ib. Since 
any analytic function of z or of z satisfies the Shroedinger equation, assume: ip = fi = 
(z — ib)°'g{z) in the upper half plane {y — 3z > 0), where g{z) is regular, and "0^/2 = 
(z — ib)°'g{z) in the lower half plane imposed, and analytically each can be extended to 
the whole plane. Obviously /i — f2 on the real line, but their normal derivatives are 
opposite and 7^ 0. (If the derivatives were zero, the functions would vanish identically.) 
The same conclusion is reached if the solution in each half plane is itself a sum of analytic 
and antianalytic parts. 
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exp(/ pdR) ~ exp(±Q;'7ri?^/2). The minus sign is excluded by the boundary 
condition at the flux sites 



2 Basics 

Let a flux of strength a be located at the origin. For an 'electron' of charge 
-e, the Aharonov-Bohm gauge potential A and the covariant derivative D in 
the standard circular gauge are given by 

iA^,Ay) = (a/2)(y/r^-x/r^) 

= {a/2)idy,-d,)\n{r), 

V X A = 2TraS'^{x), 

D = (V-iA) = (d/dx + iy/{2r^),d/dy-ix/(2r^)) (2.1) 

For notational convenience, we are taking e/i = 1 so that integer a corre- 
sponds to a quantized flux. In terms of the complex variables z — x + iy and 
z — X — iy, they become 

x={z + z)/2, y = {z- z)/2z, 
d/dz = {l/2){d/dx - id/dy), 
d/dz = {l/2){d/dx + id/dy), 
{d'^/dx' + d^/dy'') = A{d/dz){d/dz), 
Az = ia/2z, Az = —ia/2z 

D, = d, + a/i2z), D-z = dz - a/{2z) (2.2) 

The basic Schroedinger equation reads, in time dependent (energy Eigen- 
value) forms, after redeflning mE/2 = k'^/i — > E 

{{D,,D,}/2 + dr)^ = 

{{{D,,D,}/2 + E)^ = 0) (2.3) 

with the condition that ijj is one- valued and finite. (We will not necessarily 
require flniteness of derivatives, though.) Here a Euclidean time is used for 
later convenience, z and z are regarded as independent variables spanning 
a complex 2-dimensional plane. The physical space is its subspace where z 
and z are complex conjugates of each other (real x-y plane). Hereafter we 
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will refer to this situation as "on shell". Since Eq.(2.3) has independent first 
derivatives in z and ^, the actual physical space also includes its tangential 
neighborhoods. After solutions have been written down in z and 2;, however, 
we can stay on shell by reverting to x and y. The ordering of Dz and may 
be ignored since their potential noncommutativity at their singularities will 
not arise. Prom Eq.(2.2) we see that Eq.(2.3) is invariant under the operation 

2; <-> q; — > —a (2.4) 

It corresponds to a reflection y —y. The interchange of z and z is also 
effected by complex conjugation of the wave function which corresponds to 
time reversal. 

We next ehminate from Eq.(2.3) the gauge potential by a singular gauge 
transformation G: 

^ = G^, G = {z/zY''^, (2.5) 
so that -0 now satisfies a free Schroedinger equation. 

(9A- + W = 

((9A- + £^)V' = 0) (2.6) 

{A is pure gauge if its z{z) component is a function of z{z) only. See Appendix 
1.) Since the original is one-valued, 1}) must be singular in such a way as to 
cancel the singularity in the gauge function in Eq.(2.5). This means that 'ip 
must be defined as a free wave function on a Riemann surface with a branch 
cut running from to +00, and on shell there is a phase change as we go 
around the origin once: 

'0(6' — 27r) = riip{9 = 0), , r] — exp(2i7rQ;), 

^(^ = ^ = 0) = (2.7) 

and subject to conditions of finiteness at infinity on shell. In the following 
we will mainly be concerned with the eigenvalue equation. An elementary 
solution of Eq.(2.6) is 

exp(^t - Ez/t) (2.8) 

where t and i are complex momenta t. The general solution of Eq.(2.6) 
that satisfies the boundary conditions may be built up as a superposition of 
elementary solutions 

= [ ex-p{zt - z/t)f{t)dt (2.9) 
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with an appropriate choice of the function f[t) and the integration path in 
such a way as to satisfy Eqs.(2.7) on sheU. We recognize the famihar integral 
representation of Bessel functions. 

Now assume for the time being < a < 1 , and let 

/(t) = t-"— 1, n = 0,1,2,... (2.10) 

and choose the contour C of integration to be 

C : (U(?R:(zt) < 0)oo, (0+)), (2.11) 

which goes around starting from, and ending at, infinity in the direction 

U such that the integral is convergent. This means that, as the phase of z 
rotates by 2tt, the contour integration will have to make a counter-rotation 
by —271, so the factor f{t) will yield a phase factor e^*'^" as is required by 
Eq.(2.7). Furthermore, ip' stays finite on shell when z and z go to infinity, 
and also at zero since 



^-a-L-n^^ = 0, n + a>0 (2.12) 

These properties can be made more explicit by a change of integration vari- 
able, t ^ E/z 

^ = ^"+" exp{Et - zz/t)t-"-''-^dt (2.13) 

up to a constant factor, and where the contour C is now {—oo, (0+)). So 
Eq.(2.13) gives a desired set of solutions for any a and integer n = 0, 1, 2, . . . 
that satisfy n + a > 0. In the case n -|- a < 0, we can choose a new contour 
around 0: 

C" : (0 X U{m/t > 0, (0+)), (2.14) 

i. e. one that starts from the origin in the direction U, and comes back after 
encircling it clockwise. This is equivalent to the conjugate form (f-type) of 
ijj according to Eq.(2.4), which can be converted to the original contour by 
the substitution in Eq.(2.10): t E/t, resulting in a new /{t)' — 
Eq.(2.13) then becomes, up to a constant factor, 

^ = I exp{-Et + zz/t)r+''-dt, n + a<0 (2.15) 

In terms of Bessel functions, the z- and z-type solutions lead to solutions in 
the original gauge: 

n + a>0, z-tjpe : ^ = {z/z)''/'^Ja+n{2{EzzY^^), 

n + a<Q, z-type : * = (^/^)"/V|«+„|(2(£;^^) ^2) (2.16) 
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Thus the physical states are labeled by jo; + n| and an integer angular mo- 
mentum n = 0, ±1, ±2, ■ ■ ■. The angular momentum operator L is given by 
zdz — zdz- Eqs.(2.13) and (2.15) show that the nontrivial phase information 
in ip is carried respectively by z and z. Applying the operators dz and dg on 
a ip serve to change n by ^1. This may be repeated any number of times as 
long as the wave funtion remains finite, thereby generating a complete set of 
states. 

Some remarks are in order here about the problems associated with non- 
single valued wave functions and the limitations in taking derivatives. The 
first problem poses a difficulty in defining a complete set of states. Although 
scalar products are independent of the choice of (unitary) gauges, it does not 
make sense to expand a function with a branch point at z = a in terms of 
functions with a branch point at z = b. Similarly, the meaning of going to 
momentum space by Fourier transforms is not clear since it will depend on 
where a cut is made, and is diffucult to evaluate. We should, therefore, go 
back to the original non-singular gauges to build a Hilbert space of physical 
states. Wc may, however, work with non-single valued functions in enumer- 
ating the states since there is a 1-1 correspondence between the two sets. 

The second problem, on the other hand, causes a problem in defining 
higher powers of mementum as observables, e.g., the finite translation oper- 
ators in power series. For convenience, however, we may pretend to ignore 
these problems, and use such formal expressions as {ik)°' to mean — 
when applied to Eq.(2.8) in momentum space. 

3 Scattering 

Since the solutions Eqs.(2.16) are Bessel functions of non integer order, we 
expect that the scattering amplitude of an electron by the flux can be ob- 
tained by a superposition of all solutions labeled by n, which would form a 
complete set. However, since they all vanish at the origin, they form a com- 
plete set only in the punctured plane — 0, hence it would not be possible 
to construct a plane wave in with this set of functions. What is missing 
is an equivalent of the s wave component that remains flnite at the origin, 
a situation analogous to the case of hard sphere scattering. But it is also 
different because of the long range nature of the gauge potential as in the 
case of the Coulomb scattering. 
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A proper scattering amplitude can be constructed in the following way. 
Consider the integral. 

V' = / exp{Ezt - z/t)f{t)dt, 

J c'i+c'2 

fit) = {t/UnT/{t - tin): Un = lE^^^ (3.1) 

tin is the momentum of the incoming plane wave moving in the x direction: 

V'm = ex-p{iE^/^{z + z)) = ex.-p{2iE^/^x) (3.2) 

The contours Ci and C2 are defined as follows. [Fig. 1] Denote by Rq the refer- 
ence Riemann sheet, and by i?„ the sheet reached by going counter-clockwise 
around the origin n times. For a given position z = r exp{i6) = rU, Ci starts 
from 00 X f/^, goes counterclockwise around the origin and to end up at 
00 xU^ on 5*1. C2 starts from x , goes clockwise around the origin, but 
avoiding to end up at x C/^. 




Fig.l Contours for scattering 



Some properties of Eq.(3.1) are obvious: 

1) The integral converges for z = z = giving ^ = 0. 

2) When a is an integer, the two contours can be detached from their limits 
to form a closed circle around which then picks up the incoming wave 
%l]in only, so there is no scattering. 

That the formula gives the correct scattering amplitude can be shown by 
computing its asymptotic form by the saddle point method. First take the 
general formula Eq.(2.9). Its saddle points are determined by the extrema to 
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of S{t), and performing a Gaussian integration around them in the direction 
of the steepest descent. In the case of scattering, 

S' = z + Ez/tl + a/to - l/(to - tin) = 0, 

S" = -2Ez/tl - a/tl + l/(to - knf (3.3) 

For large E^/'^\z\ = E^^'^r >> 1, there are three saddle points on a circle of ra- 
dius E^/"^: one near the pole tin representing the plane wave, and a conjugate 
pair t+,t- representing radially outgoing and incoming waves respectively. 
(tin and t± do not coincide except when 9 = or n.) But the contours. Ci 
and C2 are such that the incoming component is canceled, leaving only the 
outgoing one. For the radial components, 

t± ~ E^/^exp{±ie), 

iE~^/^ exp{±2iE^^^R - iaO ± in / A) / {±e-'^ - 1) (3.4) 

The Gaussian integration around t± yields 

exp{±2iE^/^R - i{l/2 + a)9 ± iTr/A))/ R^^^ 

/sm{9/2) ioT^+, 

/cos(^/2)) foTij_ (3.5) 

These should also be multiplied respectively by an additional factor depend- 
ing on the contour: —r){r)~^) and —1(1) for Ci(C2). So when the two contri- 
butions are added, and after reverting to the original gauge, the scattering 
amplitude becomes 

ipscat ~ sin(7ra) exp{2iE^/^R - i{l/2 + a)9 - m/A)))/{R^^^ sin(^/2)) (3.6) 

which reproduces the result given by |^. The formula is valid except for a 
near forward region where the width of the Gaussian integration ~ R^^'^E^^'^ 
becomes larger than \tin — t+\. (see [] Ruij senaar s |TT| S ako da) . 

4 Presence of a magnetic field 
4.1 Basic formulas 

In this section we allow for the presence of a real magnetic field in addition 
to the fluxes. The vector potential for a constant magnetic field 7 is 

A, = 7^/2, A, = -^/z/2 (4.1) 
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which is to to be added to the Aharonov-Bohm potentiaL After gauging 
away the latter, we have 

= E^,H = -{{d, - 7^/2), {d, + iz/2)}/2, o (4.2) 

Assuming 7 > 0, make a further transformation 

\1> = exp{-jzz/2)4), H'^j = Eip, 

H = (9,7z)9,- + 7/2 (4.3) 

Since dz commutes with H, an elementary solution is then 

^p = exp{zt){t--fz)^^^-^/^, (4.4) 

where t is a complex constant. If there are no fluxes, ip must be single- valued 
and finite everywhere, i.e., the second factor must be a polynomial, which 
leads to the familiar Landau spectrum: 

E/7 - 1/2 = n = 1,2- ■ ■ (4.5) 

so that 

= exp(-7zz/2 + tz)(t-7z)" (4.6) 

For n > 0, ip vanishes at 7^ = t, and goes to zero at large \z\. The states are 
labeled by the level number n, and a continuous complex parameter t which 
is the analog of the continuous real parameter in the Landau gauge. (The 
true energy is E/4m in our convention.) Here t/^ corresponds to the center 
of the classical Larmor orbit with angular momentum —n < 0. Note that 
[t — 'jz) is a raising operator for the level number n, hence < t >=< >. 
We will conveniently refer to this factor as raising factor. 

The space of all n's and all complex t's is overcomplete on the physical 
shell, unlike in the Landau gauge. (In the circular gauge, the radius \t\ and 
n label the states 0].) But by expanding Eq.(4.4) in powers of t, or by the 
Fourier integral 

ijn,l = f i^nt-'-^dt, l>0 (4.7) 

we get a complete set centered at zero and labeled by two integers n and I, 
giving a total angular momentum j = 1 — n. Each solution contains a factor 
z^ times a polynomial of order n m z and z. Eq.(4.7) allows an intuitive 
interpretation that these states are an epicycle-like superposition of Larmor 
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orbits, with 'orbital' angular momentum /, along another circle of a certain 
radius (see below). When 7 < 0, clearly we may take the conjugate solutions 
to the above z <-> z, 7 — — 7,j — > —j, and the energy levels are given in 
general by -E/I7I = n + 1/2. It is also instructive to examine the limit 7 — *• 0, 
where we should recover the original free field results. Since the energy 
spectrum is proportional to 7, however, in general it will be driven to zero in 
this limit, as may be seen from the elementary solution Eq.(4.6) which does 
not have a limit E = n/'y except for n = 0, E ^ 0. But supplying a factor 
i.e. dividing the raising factor by t, we get 

lim exp{--fzz/2 + zt){l - 72/^)^/^"^/^ = exp{zt - zE/t) (4.8) 

which reproduces the free field exponent. 



4.2 A flux in a magnetic field 

There are two ways to introduce a flux of strength a. One is to change n 
to n + a in Eq.(4.6) (provided that it is > 0), and to let it be the gauge 
transformed ip. It represents a flux located at t/7, and the energy is shifted: 

E/7- 1/2 = n + a > (4.9) 
By expanding Eq.(4.9) in powers of t, or taking the moments 

i^n,l = j ^Pn+a,tt-'-^dt, < I < k = [u + a] (4.10) 

we get Iq states with the same shifted energy. The restriction / < Iq arises 
from the condition ijj{0) = 0. This is a z type representation. Q Both Eq.(4.9) 
and Eq.(4.10) are now nonpolynomial functions in z. 
For I > Iq, the proper formula is given by 

Jc 

E/7-l/2 = n, C = (oo,(0+,7^+)) (4.11) 

■^It can be turned into a z-type representation if we change the sign of 7, and let the 
path be C = (00, (0+)). The energy ^ E/^ remains positive. See section 6.3 
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which is a -2-type representation. The energy does not get shifted for them. 
For n = 0, in particular, if) is the same as the zero energy hmit of the field- 
free product of ZiS, apart from an overall Gaussian factor. The 
reason why a nonpolynomial raising factor in the integrand is not allowed 
is that, at large distances, ip becomes ~ exp(+7z^/2), as may be seen by 
finding the saddle points. It is not possible to avoid the exponential blowup 
by changing the sign of 7 or energy. The problem is avoided only if the 
factor is a polynomial. (According to the previous subsection, it can also be 
seen that only unshifted levels have zero field limits.) Obviously these same 
results will hold for 7 < 0, for which the conjugate representations are used. 

Summarizing, a number Q < Iq = [n + a] oi states at level n have their 
energy (more precisely -E/I7I — 1/2) shifted to n + a > 0, and the rest is 
uninfluenced by the flux. Consider in particular the ground states n = 
with \al\ < 1. If a > 0, 7 > (or if a < 0, 7 < 0), i. e., the magnetic field 
and the flux are parallel, then the energy of one state is shifted upwards by 
a. If, on the other hand, the field and the flux are antiparallel, none of the 
ground states can change energy; the energy can only go up, but not down 
below the zero point minimum. 

The following semiclassical argument explains the difference between shifted 
and unshifted levels (see also [|14|)- The state 'iljn,t around t occupies a disk 
of radius R ~ (72/7)^/^. The orbital radius of the 'epicyle' state il)n,i is 
r ~ IjE^I'^ ~ l/{'ynY^'^. So the latter state will cover the flux at the ori- 
gin only if i? < r. As for the asymmetry in the ground state, write the 
Hamiltonian as 

H = p^ -L{a/r'^ + -i) + {a/r + -irf (4.12) 
For L = 0, the minimum of the last term is 07 if 07 > 0, and if 0:7 < 0. 

5 The multi-flux problem 

5.1 Multiple integral representation 

At zero energy, the general solution of the Schroedinger equation is either 
analytic or antianalytic. The one vertex solution reduces to ^/^ ~ or z". 
More than one vertices can be easily accommodated by taking products of 
n such 2;-type or ^type factors, but we cannot form general mixed products. 
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When E > 0, the solutions will involve both z and z, but by continuity 
this property will persist. In the following we assume nonzero energy. Our 
goal is to find a solution which has proper monodromy properties around 
each flux:Mjt/' = rjiijj, where Mj denotes the monodromy operation around 
flux i, giving a phase rji — exp(i^j). Since, however, the radial second order 
differential equation obtained after removing G has singularities only at 
and oo, it is not clear how to introduce more flux singularities. 

So let us flrst try a peturbative solution in energy for two fluxes. Denote 
their location and strength by (61, ai), (62, ^2), and write z — bi = zi,z — b2 = 
Z2 for short: -0 = E ^^""^z , z) , = (^2)"'- Then 

£;0('^+i) = -a^a^-V;(") (5.1) 

ijj^^^ can be obtained by integrating ip^^^ over z from 61 (&2), and over z 
from bi (&i). This -0*^^^ vanishes at both sites, but it does not have a proper 
monodromy at the second (flrst) site. We can make behave correctly 
at both sites by taking ij/'^^ = {zi)°'^'^^{z2)^'^^^'^z with integer rii > 1. ij/*^^ 
is then obtained by differentiating But to generate a series to the A^th 
order, we have to start from an A^th order term of sufficiently high power 
Hi > N, and we lose copntrol of the perturbation expansion. 

Next we look for a nonperturtive solution in terms of a double integral 
representation for two fiuxes. The gauge potential and gauge function are 
respectively a sum and a product of two components: 

= ai/zi + a2/z2, = A^, 
G{z) = {zi/z,r/\z2/z2r/' (5.2) 

Consider a solution of the Schroedinger equation of the form 

= exp[2;iti + 2:2^2 - E^/'^{ziSi + ^2S2)]/(^i, ^2, Si, ^2), 
{ti + t2){si + S2)^E (5.3) 

The idea is to integrate over one of ti, si and one of t2, S2 independently to 
satisfy proper monodromy, but the constraint precludes mixed types involv- 
ing t and s, and it will be sufficient to examine the z type only. Since the 
third variable turns out superffuous, we are led to the integrand 

= exp[zih + Z2t2 -{z- b)E/{h + t2)]tT^t^'''-\ (5.4) 
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with an arbitrary b. The integration contours (Ci,C2) are formally taken 
respectively around to U{—6i)) and U{—92)- It is in general not possible, 
however, to avoid the vanishing of ti + 12 by a proper choice of the contours, 
and this would invalidate the formula. As Ci or C2 makes a 27r rotation, 
there always occur crossings of — Ci and C2, and we have to set up a con- 
vention for one of the t's to make a detour around the other. A change of 
the crossing point leads to a change in the integral, and after making a 27r 
rotation, it will cause an extra change in ip beyond a phase factor. A de- 
tailed study shows that the extra piece is a regular function of z and spoils 
the monodromy properties. Therefore it is not possible to satisfy the two 
monodromy conditions simultaneously. Later an alternative form of single 
integral representation for dynamical fluxes will be found to solve the prob- 
lem. But we will first show next that the multiple integral representation can 
be used without difficulty if a magnetic field is present, athough the solutions 
obtained are of the type that do not have free field limits. 

5.2 Solutions in a magnetic field 

Unlike the field-free case, it is easy to generafize the -z-type representation for 
unshifted levels, Eq.(4.11), to N fiuxes with the multiple integral method. 
The solution is given by 

= exp(-7^^/2) / ■ ■ ■ exp(5: z,U){T - ^zf nr"^~""'^^0, 
T = = + > 0, i = \...N (5.5) 

The factor (^^(-Zjti — 7^)" can become zero, but this does not cause the earher 
problem since it is a polynomial in the tj's containing powers up to iV. The 
condition on the a^'s insures that the overall power of each ti is < 1. For the 
ground state n = we have 

i^{o.'}fi^V^zt' l{ol)\ *o = exp(-7^^/2)V'o (5.6) 

The higher levels are generated by applying the raising operator dz — z'y n 
times to '?/^{a'+n},o- Thus for this class of solutions the Landau levels are not 
affected by the fluxes. The wave function is a homogeneous polynomial of 
order n made up of z^ and the Zi%. 
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5.3 Zero field limit 



Earlier we found that the multiple integral representation failed in the field- 
free case, so it is an interesting problem to study the limit of letting 7 go 

to zero, 7 = E/{n + 1/2) — > 0, keeping energy fixed, or n — >■ 00. We will 
evaluate the formula Eq.(5.5) for N fluxes in the n ^ 00 limit by finding the 
saddle points with respect to each ti. 

U = (a[ + n)/(zi + n/(T - Ez/n)) T,T ^Y.^, ^ NT (5.7) 

which is a contradiction for > 1. A more careful analysis shows ~ 0{n). 
The wave function is then 0(n^^~^^"), and it is driven to zero. For N — 1, 
on the other hand, a consistent limit exists as was observed in section 4. By 
writing 

{a[ + n)/ti = {zi + n/ti + Ez/tl),ti = a[/{zi + Ez/tl) (5.8) 

we get the same equation derived from the zero field solution Eqs.(2.8),(2.9). 
The origin of the problem is traced back to the fact that each factor in 
Eq.(4.16) had a power —a' — n in order to give negative powers of ti for all 
the terms ~ j^rn^^^-^n-m -j^ ^rp _ r^^Y^ and this contributed an overall power 
of Nn. 



6 Dynamical fluxes 

6.1 Wave functions for the charge-flux system 

The difficulties of the general iV-fiux problem can be resolved if we treat the 
fluxes themselves as dynamical objects with their own kinetic energies, for 
reasons which are not immediately clear. First we note that, when the fluxes 
are made dynamical objects, they are under the influence of a magnetic (dual) 
vector potentials produced by the charges which, however, are of the same 
form as the electric vector potentials produced by the fluxes and acting on 
the charges. This is obvious because, for the wave function ^0 as a function of 
the coordinates of a charge and a flux, a rotation of the charge around a flux 
should be equivalent to a rotation of the flux around the charge. (A general 
description of charges and fluxes as dual objects is given in the Appendix 2.) 
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Thus the gauge transformation G of Eq.(2.5) will remove the potentials from 
both the charge and the flux, and we will end up with a free Hamiltonian 

H = H, + H, = -d,d, - (6.1) 

where the flux is represented by coordinates 6, h and mass yU (the mass of the 
charge = 1). This generalizes to any number of fluxes as well as charges if we 
consider only charge-flux interactions, and ignore charge-charge and flux-flux 
interactions. 

The center of mass commutes with the free Hamiltonian, hence it is a 
constant of motion even with the boundary condition imposed by the flux 
singuralities since they depend on relative coordinates. So the wave function 
can be separated into that for a free center of mass motion and that for the 
relative motion, and we may deal with only the relative part. We will not do 
this here because it will complicate the formulation. 

Suppose there are N fluxes with masses /ij, i = l..A^. We relabel the 2N + 
2 coordinates of the charge and the dynamical fluxes z, ^, {fej, hi} uniformly 
as {ui, Mj}, i = 0..N and similarly the masses as {/Xj}, /io = 1- The fluxes will 
eventually be set to their classical positions 6° = u^,i ^ 0. The coordinates 
relative to them are denoted with primes: u[ = Ui — u^,i > 0). The total 
Hamiltonian if is the sum of N individual ones 

Hi = = H (6.2) 

Further define 

li = exp{-fXiUiu[/t)/t),I = Y[li (6.3) 
The J's are the heat kernels satisfying the relations 

HJi = {-^i,UiU^/e + l/t)h = -dth, HI = -dtl, {t ^ 0) (6.4) 

Here t may be regarded either as an imaginary time in a time- dependent 
Scroedinger equation, or simply as a parameter to be integrated over. ^ 

■^With the first interpretation, there are time-dependent solutions for non-dynamical 
fluxes of the following form, using only the kernel: Io{z) : tfj — Io(j){z/t), (j) — Di^iAf (^/* ~ 
ivi)°'' . The flux coordinates bi — v-it are moving with velocities Vi. But since the hiS axe 
fixed, it does not correspond to physically moving fluxes. (However, this form is of the 
same type as the function (j) given below.) Solutions with dynamical fluxes of the simple 
form exp(iPX -I- iPX) Y\{z — 6i)"S where P and X refer to the center of mass, are not 
admissible since they are not bounded at infinity. 
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The ansatz for ip is now, 

11)== [ exp{Et)I(l){ui,t)dt, (6.5) 
Jc 

with some contour C. Now apply H to it, making use of Eq.(6.3): 

Hij = J exp{Et)[{-dtI)(j)-Y.^/f^^d^'Jdu..<l>]dt 
= + J {- exp{Et)[dtI)(l) + lY.iu^du,(l)]dt 
= E^ + J exp{Et)DI(j)dt, 

D = tdt + Y.^idu^ (6-6) 

So ip will be a solution provided that Dcp = 0, i. e. 4>{ui, t) is a homogeneous 
function of degree zero. We will choose 

<P= U {{uo-u,)/t)< (6.7) 

i=l...N 



with a contour C = (— oo, (0+)). Thus 

= Jexp{Et-J2u,u'Jt)l[{{z-b,)/t)<dt/t' 

R = {J2fiiU,u'^'^^ (6.8) 

up to a numerical factor. Its conjugate representation is obtained by replacing 
the Ui,a' by Ui, —a'^, and changing the contour to (0, (0+)) (or equivalently 
t — >■ 1/t with the original contour). The above ip obviously has the right 
monodromy properties. When the flux masses become large, /ij — oo, the 
Bessel function will become small and oscillate rapidly except at their in- 
tended 'classical' positions where E^^'^fii >> \u[\ {i ^ 0), R^ z{z). For 
large values of \z\,il) behaves like a single- flux solution with strength Yl, a'+N . 
Q In momentum space, if) has a simple but symbolic form 

^(fc„fc,) = ^(^ hh/^Ji, - E) \{{i{k, - ki))< (6.9) 

^There is a freedom of including an optional factor z/t representing a nondynamical 
flux at the origin. When no other fluxes are present, this is equivalent to Eq.(2.13). We 
also note that even a single flux, not located at the origin of coordinates, does not have a 
static representation. 

'^For \z\ » \bi\, (f) satisfies the truncated form of Eq.(6.6): D'cf) = {zdz + tdt)4> = at 
the saddle points, which assures, as may be seen easily, that -0 is a semiclassical (Hamilton- 
Jacobi type) solution of H^- 
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The Bessel function in Eq.(6.8) corresponds to a solution of the wave 
equation in a (2A^ + 2)-dimensional space in the so-called hyperspherical har- 
monics expansion. In the absence of the factor 0, the general solutions in 
this space are spanned by Bessel funcions times the (2A^ + 2)-dimensional 
spherical harmonics. They can be generated by applying du^ and which 
commute with one another and with if, repeatedly to the s-wave solution 
Jjy/R^. Obviously each operation respectively generates a negative or pos- 
itive angular momentum around the origin in the 2-dimensional subspaces 
spanned by Ui,Ui. In the integral representation of Eq.(6.8) without 0, ap- 
phcation of du^ (dui) will bring down a factor Ui/t iui/t). Application of a 
monomial of these operators will generate a polynomial in the m's and m's, 
and this way a complete set of states labeled by the 'magnetic' quantum 
numbers of 0{2N + 2) and a Casimir invariant (the degree of the monomial) 
will be generated. Fractional factors in are of a similar nature, but the 
power of angular momentum around each flux are correlated to be all posi- 
tive for tp, or all negative for its conjugate representation. Nevertheless, by 
applying these same operators to the of Eq.(6.8), new admissible states 
can be created as long as the lowest power of each factor {z — hi) is positive. 
Together with their conjugates, these will presumably span a complete set of 
states. 

The solutions of the type given by Eq.(6.8), however, are not the natural 
ones when considered as wave functions for each of the individual particles. 
Because of the factor l/t^"*"^ characteristic of the (2iV-|-2)-dimensional space, 
they go like l/R^^"^^"^ ~ l/\ui\^^^/'^ when one of the Wj's becomes large, 
contrary to the expectation that they should go like in the physical 

2D space. The extra power of 1/t which is responsible for this is a quantum 
effect associated with the dyamical degrees of freedom of each flux. If a 
flux is to be pinned to a fixed position, its degrees of freedom will have 
to be frozen. By the same token, if two fluxes are to be merged, those of 
relative motion will have to be suppressed. But the flux velocities d^/ ix 
(unless evaluated at the origin) are independent of /i, which was necessary 
for construction of the factor 0. So the derivatives, i.e. the dynamics of the 
fluxes cannot be ignored at the level of the Schroedinger equation even for 
large masses. Herein lies the need for treating the fluxes as dynamical, and 
the difficulties of pinning them down. 

^However, in the opposite limit of the single charge becoming infinitely heavy, the wave 
function of the fluxes would reduce to a product of the individual ones as in the case of 



18 



A sensible way to recover the proper 2D asymptotic behavior may be 
to produce highly polarized states in the direction of the particular pair 
of coordinates Ui,Ui, i. e., those consisting of high powers of homogenious 
monomials Ui/t and Ui/t. But this cannot be a finite series since they increase 
powers ofl/t rather than reducing them. Consider then applying the Green's 
function 

P = 1/{E, - H,) = l/{d,d, + E,), E,<E (6.10) 

to the ip of Eq.(6.8). P commutes with H so Pip is again a solution, and 
should go like cos{2Ey^\z\ + 5)/\z\^^'^. Since the source ip is concentrated near 
the origin (~ l/l^l^"*"^/^), it is clear without explicit computation that Pip 
will go as 1/1 2; I In order for this to be an exact solution for all z, however, 
P must be a Green's function in which the flux singularities are already 
incorporated. Such a Green's function can in principle be constructed from 
a complete set of states, but it is not practical to carry it out. 

Instead, we can define an effective Hamiltonian for which solutions with 
the desired porperties can be written down easily: 

N N 

= -z~^d^ Y Uidui - z'^d^ Y ^i^ui 
1 

N N 

= -z~'^d^Y^i^y-i - ^'^^^IZ^i^y^i (6-11) 
1 

with the stipulation that elementary solutions depend only on either the 6's 
or the 6's, and satisfy proper monodromy and boundary conditions. Indeed, 

^ = J(-zz/t + Et)Y[{z-bi)/t)<dt/t 

= Jj^^,{2Ey'\z\)l[{z-k)/\z\)< (6.12) 

and its conjugate form clearly satisfy Eq.(6.11), and have the right aymptotic 
behavior. Being linear in the derivatives, we may say that the flux coordinates 
retain only partially their dynamical status. The general properties of Hes 
will not be exmanined here, but Eq.(6.12) will be adopted in subsequent 
sections. 



many charges in the presence of a single static flux. 
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6.2 Solutions in a magnetic field 

The Hamiltonian for a charge in a magnetic field is 

H, = -{d,-z-i)d, (6.13) 
and a solution which has a zero- field limit is, according to Eqs.(4.8), 

tl) = j ex^{zt){l-z-f/tfdt/t 

= J exp{-zz/t){l + -itydt/t = J IJt/t 

h = exp(-^^/t)(l + 7t)7i (6.14) 

The last expression was obtained by a change of variable t — > —z/t. It 
satisfies the relation 

{H, - 717)7, = -dt{{l + -it)Q (6.15) 
By substituting the new J,, the analog of Eqs. (6.4-5) becomes 

if) = j I(pdt, 

(H-E)^ = J -{d,{l + ^t)I)<f> + J2{^,/t)du,<l>, 

ID'(j)dt/t, D' = {1+ ^t)tdt + Y,Uidu, (6.16) 



/ 



The solution for (f) is then obtained by substituting t/(l + 7t) for t in 
Eq.(6.5), and becomes (dropping exp(£^i) in favor of (1 -|- 7^)") 

j exp(-^x.,iZ;A)(l/i + 7))^(l + 7^r 

\{{{z-hi){l/t + ^l))<dt/t (6.17) 

The integrand has a branch point at t = of order — — ^a' and one at 
t — — 1/7 of order N + n + YjOi' . We can adopt a closed contour (— I/7-I-, 0-I-) 
around them. When the integral is evaluated by expanding the exponential, 
we get a series of Beta functions 

{jrYB{-N - ;^a'-m,Ar + n + ;^a' + l)/m! (6.18) 

which is zero if m > n + 1, so it is a polynomial of order n. In the limit 
7 — > 0, n = £'/7 — > 00 for fixed E, (1 -|- 7^)" becomes exp(£'t), the contour 
becomes (—00, (0-|-)), and Eq.(6.17) reduces to that for the field-free case. 
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6.3 Shifted Landau levels 



In section 4.2 we showed that there are shifted Landau levels for a single 
static flux. Eq.(6.17) serves a basis for incorporating N dynamical fluxes. 

Assume this time that 7 < 0, replace n in Eq.(6.16) by —n — J2 a' — 1 so that 
the energy E/'-f = (n + + l/2)/|7| is still positive. The singularity at 
— 1/7 becomes a pole of order N — n. With a change of variable t = s — I/7 
the full wave function ^ can be turned into 

^ = exp{—'~fzz/2)ijj 

= / exp(-r'V(s- 1/7)- 7^^/2)7^-"-' n((^-^0"^ 

(s - l/7)-^-S"'-^s^-"-^ds (6.19) 

which is nonzero if n > iV, or if n > after N is eliminated. \1' converges at 
00 since 7 < and r' — zz effectively. The energy is n + a'^. This, then, is 
the generalization of the shifted Landau levels, Eqs.(4.9), (4.10). For a single 
static flux plus N dynamical fluxes, the energy gets shifted by all the flux 
strengths simultaneously. We do no know if partially shifted states by some 
of the a'^s can be found. When 7 — > 0, n — > 00, the integral is pushed to 00 
and disappears, meaning that the shifted states have no free field limits. 

7 Large N limit and flux medium 

According to the results of Section 6.1, we will effectively regard the fluxes as 
fixed objects and deal with wave functions of the type given in Eq.(6.12). The 
order of the Bessel function representing ijj then goes up with the number 
of fluxes as with all a's being of the same sign > and adding up. 

Suppose the fluxes of equal strength a are distributed with a density l/a^ up 
to a large radius Rq. The total flux is an{R/ay. The angular momentum L 
of the wave function at radius R < Rq {R > Rq) is aTi{R/aY {a7r{Ro/ay). 
Its azimuthal momentum is L/R = anR/a? {aTr{Ro/a)'^ /R), so the radial 
momentum for large R : airR^E >> 1, is ~ ^iamRj , and the radial 
"phase" is / i'PRdr ^ ±aTi{R/af /2 {±{Rq / af \q. R) . For R < Rq, the wave 
function will behave like Gaussian or anti-Gaussian. Below we will investigate 
it in more detail. 

Assume that all fluxes have the same strength a, and form a regular 
square lattice of unit size at sites bi,i — m + in, although these assumptions 
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are not essential. The gauge potential A^, and the gauge function G are 

A, = (a/2)^l/(z-6,), 
A, = -{a/2)J2l/{z-k), 

G = m^-'bi)/i^-^i)T'" (7-1) 

As the sum goes to infinity, they fail to converge, hence an extra pure gauge 
subtraction terms become necessary (except for i — (0,0)). The modified 
forms are then 

A, = (a/2)C(;.),^-=(-a/2)C(^), 

G = [a{z)/a{z)r'\ 
a{z) = zX{{{l - z/hi)eMz/h + z^/2hD) (7.2) 

Here C and a — exp(/ C,dz) and are the elliptic functions of Weierstrass with 
full periods 1, i. Because of the subtractions, they are not periodic but satisfy 

(^{z + m + in) = ({z) + n^m — in), 

a{z + m + in) = a{z)exp[Ti{m'^ + n'^)/2 + iTi'^mn/4: + zTi{m — in)ll7.3) 

The leading asymptotic factors are insensitive to the regular lattice structure, 
and may be derived by a continuum approximation. If the flux density is 
l/a^ 

'^l/{z — bi) J l/{z — r exp{i9)rdrd6/a'^ 

= TT\z\'^/za'^ = nz/a^ (7.4) 

The contributions to the integral over b from r > \z\ vanish because the 
phases of the contributions from points at radius R > \z\ from the origin 
cancel. Integrating this from to 2; along the radius, we get 

^n(^ - b^)/b^) =(^R-- exp(7r(|^|/a)V2) (7.5) 

These non-analytic expressions are a result of confusing z and the cell coor- 
dinates (m, n)a. We can also use this method to deflne a cutoff version of C, 
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and cr: 



(fi^{z/a) ~ TTz/a if \z\ < Rq, 

= if \z\ > i?o; 

(TR^{z/a) ~ exp(A^/2) if \z\ < Rq, 

aR^{z/a) ~ exp(iV/2) = vri^^if \z\ < R^, 

~ exp(A^/2)(z/i?o)^ if \z\ > Rq, (7.6) 

We can now substitute them in Eq.(6.12), taking the origin of coordinates 
at the center of the lattice. Q Then = cxij/t"^. First, the value of ip at the 
origin can be exactly calculated, and yields ip{0) = 1/N\. For the asymptotic 
behavior, we can use the the saddle points to get 

ij ~ C{to) exp{-E\z\yto + h)^'"'an,{z/ar (7.7) 

where C(to) is the contribution from the Gaussian integral around t^, and 

to = ±iE^/'^\z\ if N»E\z\^, 

iIj ~ eM^'^^E^'^\A+Na/2){z/ ±iaE^'^\z\)^''/{±iE^I^\z\f^] 
to = Na if ERl < E\z\'^ < N, 

~ A^i/2gxp(3Ara/2)(z/aA^a)^" 

= exp(Ara/2)(z/a)^"/Ar! 

= exp(A^a/2)(z/a)^X0); 
to = Na if 1 < E\z\'^ < ERl, 

~ exp(a7r|2|V2)(z/|2|)""l"lV(0) (7.8) 

The first form is the true asymptotic limit. The second is valid in the near 
outside region \z\ > Rq. The third applies to the inside of the flux medium. 
These are all consistent with the physical argument given at the beginning 
of the section. |^ From above we see that the wave function inside the flux 
area grows with the distance in an anti-Gaussian way like exp(a7r|2;p) inside 

possible alternative choice for <j) is cj) = aR„ (z/a) joji^ (It) with some / and the contour 
running outside of the fluxes: It > Rq. But this leads to the same asymptotic results as 
those given below, up to a normalization. 

®In the exterior region, cj) satisfies D'(f) = zdzcj) = —tdtcj) = aN, so ■0 is a semiclassical 
solution of Hz (footnote 5), whereas D'(f> is not zero in the interior region. (Of course the 
full D(f) = is always guaranteed). 
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the medium, then with a power hke to distances E^/^\z\ > Na where 

the true asymptotic behavior sets in. In other words, the wave function is 
pushed out ot the flux region and beyond (unless the energy is so large as to 
skip the intermediate regimes). 

8 Summary and discussion 

In terms of physical concepts, the Aharonov-Bohm problem for many fluxes 
is a straightforward extention of the single flux case, and is not expected to 
contain anything conceptually new. Yet in an attempt to construct explicit 
solutions, we have found some unexpected features. In the idealized A-B 
problem of an infinitely thin flux, the wave function of a charge in its vicinity 
is described by Bessel functions whose order is shifted by the flux strength n 
to |q; + n|. Since the latter has to vanish at the flux, the component which 
would be an s-wave in free space is pushed out irrespective of the sign of a. 
This effect will magnify in a medium made up of fluxes. 

To handle the general case of many fluxes, the gauge potential may be 
removed by a singular gauge transformation. The Hamiltonian is reduced to 
that of a free particle except for boundary conditions. Around each flux i the 
transformed wave function must develop a branch point in one of the complex 
coordinates, with a positive fractional power so that it has the form 2;"'+"% or 
^-ai+rii^ times a holonomic function in z and z. In particular, at zero energy 
the general solution is an arbitrary sum of a function of z and a function 
of z, and each part must carry all the branch points. In other words, the 
fractional angular momentum around each flux, as opposed to the intrinsic 
flux strengths must all be either positive for the analytic, and negative 
for the anti-analytic solution. At nonzero energies, the wave function of 
eithertype contains both variables, but by continuity the above qualitative 
features do not change. (This turns out to be true even in a magnetic fleld 
where the spectrum is discrete.) The total angular momentum of a wave 
function around a large circle containing a cluster of fluxes is then a sum 
of positive or negative numbers which will grow with the number of fluxes, 
resulting in the expulsion of a charge from the interior of the cluster. This 
is because the wave function has to vanish at each flux site, whether the 
fluxes add up to an integer or not. Only if an individul flux is integral, will 
it not contribute to the effect. This conclusion does not require an explicit 
construction of solutions. 
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These general results do not seem to be due to our reliance on the singular 
gauge and analyticity, or the simplifying assumption of infinitely thin fluxes. 
To see this, first consider a flux ai at the origin in the original gauge so a 
solution may be expanded in Bessel functions J\ai+n\ of positive fractional 
order. If there is a second flux 02 at point P, they have to be re-expanded 
in Bessel functions J\a2+m\ around P. But the complete set of eigenfunctions 
for the latter covers only the space — P. Therefore we have to start from 
the beginning with a function that has the proper behavior at both sites. 
(The general addition formulas for Bessel functions show that all positive and 
negative powers and with the same a are involved in the re-expansion around 
P if a is not an integer. The integer case is special since J_„ = (— 1)"J„.) 

It is also instructive to compare the A-B problem with that of hard core 
potentials. Here the similarity is that the eigenfunctions do not form a com- 
plete set in R^. The difference is that in the latter the eigenfunctions are 
Bessel functions of integer order, the s-wave included. If the radius of the 
potential is finite, the solutions involve both regular functions J„ and singu- 
lar ones Nn- This is equivalent to having both positive and negative orders 
for noninteger n, or the Hankel functions -ffo+n. In a realistic A-B problem 



where the fluxes have a finite size[|l^]|]T2|, the solution around a flux will 



be a superposition of functions Ja+n with all n, or, after gauging away the 
potential, it will have all powers 2;"+"- or This, however, would not 

alter the difficulties with many fluxes since the quasi-free wave function has 
to satisfy not only the condition that it vanishes at the boundaries but also 
that it has the right fractional angular momenta in the neighborhood of the 
fluxes, in contrast to the hard core case where only the boundary condition 
matters. 

The explicit construction of wave functions for many fluxes has also 
turned out to be possible, with or without a magnetic field. To construct a 
solution in free space, however, the fluxes have to be treated as dynamical ob- 
jects with their own Hamiltonians. Technically we may say that it is because 
kinematic correlations are induced between fluxes indirecty through those 
between the charge and the fluxes, and are independent of the masses. But 
the role of the magnetic field in this context is unclear. These are puzzling 
features that remain to be understood in simple physical terms. Explicit 
solutions with proper asymptotic behavior can be witten down for an effec- 
tive Hamiltonian which partially retains the dynamical nature of the fluxes. 
With a magnetic field, it is noteworthy that some of the Landau states get 
the energy shifted by the fluxes in proportion to the flux strength while others 
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remain unaffected. The effect is asymmetric for tlie lowest level with respect 
to the relative sign of the field and the flux. This nonlinear behavior may 
lead to some interesting consequences. 
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Appendices 

A.l Complex metric 

Let a covariant and a contravariant vector be denoted in general by v, and 
V respectively, and let the contravariant vector v = (x, y) be mapped to 
(z, z) — {x + iy, x — iy). The metric for the complex pair is then off-diagonal: 



^,fe = (l/2)[[0,l],[l,0]], g 



ik 



2[[0,1],[1,0]], 



(A.l) 



The rule for forming a vector A in the complex metric is 



A = {A', A') = (A^ + lAy, A"" ~ lAy), 
A, = {A„ A, = {A^ - iA\ A" + iA') /2 



(A.2) 



Since the gradient operator is 



grad — {d/dz, d/dz), 



(A.3) 



and the Laplace operator becomes 



= 2{d/dz, d/dz} 



(A.4) 



Rotation R by an angle 6 is represented by a diagonal matrix 



R = iv ^,V), R = {V,V V = exp(i^) 



(A.5) 
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so the outer product of two covariant or two contravariant vectors is given 

by 



V X w = -i{v''u^, -v'^u'^) (A. 6) 

It follows also that the divergence and curl (= Fzz) of A are given respec- 
tively by 

divA={dzAz + d,A,)/2, 

curl A. = -curM- = i{dzA, - d-zA,)/2 (A.7) 
The Gauss' and Stokes' formulas read (on shell) 

Vnds ~ ~^ j i^zdz — Vzdz)^ 

Vgds = f [vzdz + Vzdz) (A. 8) 



respectively. 

The Laplace operator applied to the Green's function 

G = In |z| = (In z + In z)/2 (A.9) 

is zero, but the delta function singularity at the origin of the on-shell space 
can be seen from the Gauss' formula 

j grad Gnds = (-i/2) j [dzjz - dz/z) = 2tx (A. 10) 



A. 2 Field theory description of charges and fluxes 

The fluxes can be regarded as dynamic magnetic particles. They interact 
with the charges through the Aharonov-Bohm effect, but not directly among 
themselves. A flux $ is a concentration $ = Ba? of the magnetic field density 
B over an infinitesimal area a^, its energy (per unit z direction) is a part of 
the energy density of the magnetic field plus an amount necessary to 

keep it concentrated. [Q We will assume uniformity and ignore the dynamics 
in the z direction. When a flux moves, the magnetic field inside will acquire 
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electric field components of order v/c, so {B^, —Ey, E^.) will be associated 
with a moving flux. The vector potential associated with E will be present 
in the medium, even though it may be a shielded potential. 

Let A be the usual vector potential, G and F two vector fields, J and K 
the electric and magnetic currents in 2 + 1 dimensions, so they have mass 
dimension 2. J is then not the usual charge density since the charge is not 
extended in the z direction, but may be thought of as a kink soliton d^cr of 
an extended field a. This is the only place where the z direction enters. Let 
the Lagrangian L be 

L = -J^Gi ■ Fi/2 -iG ■ cm\A + eA- J + gG ■ K + Lmat, 
Gi = D2, G2 = -D,, Gs = Hs, F^AG, 
A = diagonal matrix [A, A, (A-H) 

Lmat refers to the kinetic part of 'matter' (charge and flux) Lagrangian. Eu- 
clidean metric is used for simplicity, with 0:3 = ixo, etc. The electric and 
magnetic charges e and g are dimensionless; G and F have mass dimensions 
1 and 2 respectively. B, H, E, D have the usual meaning in the Maxwell the- 
ory except for dimensionality: the parameters k, A are respectively the analog 
of magnetic permeability and inverse dielectric constant with the dimensions 
of mass. Eq.(All) shows that the magnetic "vector potential " is invari- 
ant under the usual gauge transformation of the electric vector potential A. 
Varying G and A as independent fields, we get the field equations 

F + icuT\A^ gK, icmlG^eJ, (A.12) 

If the magnetic current K is absent, Eq.(Al2) takes the usual Maxwell form. 
If, on the other hand, k and A are set to zero, the equations show a symmetry 
under duality between J, A and K, G: 

i curl A = gK, i curl G ^ eJ, (A. 13) 

and an invariance under both an electric and a magnetic gauge transforma- 
tion: 

A^ A + grad 0, G ^G + grad x (A.14) 

(Boundary contributions are assumed to vanish uder these gauge transfor- 
mations.) The left-hand side of Eq.(Al3) are the Coulomb- Lorentz forces 
acting respectively on charge and flux, which are seen to come directly from 
the currents of their opposite numbers. 
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For static electric and magnetic sources located respectively at Xj,yj and 
Xk^Vk: Eq.(Al3) yields 

Ai, A2 = gi-v2, vi)/ {n\vf),v = {x- Xk,y- Vk) 
Gi, G2 = e{-U2, Ml)/ (7r|Mp), u= {x - Xj,y - yj) 
A+ = —iAi + A2 = g/z, A" = iAx A2^ g/z,z^vi^- iv2, 
G+ = iGi — G2 — e/w-i G- = iGi — G2 — —e/w, w — ui + iu2 (A. 15) 

Hence the vector fields eA and gG acting on charge and on flux are the same, 
and proportional to eg, and the same form of gauge transformation removes 
the potential from the charge and the flux Hamiltonian. 

If A 7^ 0, there is no duality symmetry. Conservation of the magnetic 
charge requires 

divAG' = (A. 16) 

which is the same as in the Maxwell case. The magnetic gauge invariance 
still holds if the x iii Eq. (A14) satisfies 

divAgradx = (A. 17) 

By taking the Coulomb gauge diAi + ^2^2 = 0, Eq.(Al3) leads to 

V^As = -igX-\d xK)u- AeJg, 
A-^V^i + K-'dA = -ig{X-^V2Ks + ^-'8^X2) - ieJi, 
X-^V^A2 - K-^d3Fs2 = -ig{X-^ViK3 - K'^dsKi) - ieJ2 (A.18) 

If static, ds = 0, or if k » A, the contribution to A from K, and one to G 
from J become 

V% = -XeJs-tgX-\dxK),2, 
V% = -tgd2K,,V^A2^igdiKs, 
G3 = iK^^g{d X A)i2, 

Gi = iA-ia2^3, ^2 = -iA-^ai^a (A.19) 

Since the Green's kernel is ln(r)/27r, r — ^{x"^ -\- y^), the potentials ^41^2 and 
Gi_2 are then found to be the same as Eq.(A13): 

eAa = ieg ln(r)/7r, eA,2 = -^Gi,2 = eg{d2, -di) ln(r)/7r (A.20) 
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Eq.(Al8) shows that the magnetic 'Coulomb' potential G3 acting on mag- 
netic charge qKq = —igK^ is, apart from self interaction, indeed the usual 
magnetic field G12 = = B^/n, where k is interpreted as an effective 
magnetic permeability. The term G12-F12/2 = is the magnetic field 

energy. When an external magnetic field Hex is imposed, we should change it 
to K{H+Hexf /2, which means that the flux will feel a magnetic Coulomb po- 
tential Gex- By the same token, if an external electric field is present without 
suffering shielding, the flux will feel an electric Lorentz force. On the other 
hand, the potential acting on a charge at (x, y) and generated by another 
charge at is — A(ln(r)/27r)eJ3. Since this should be the 3-D Coulomb po- 
tential l/(47rr), we must interpret A to be ~ dz-, i.e. the derivative in the 
hidden z direction at 2 = 0. f\ 
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